Momentum distribution dynamics of a Tonks- Girardeau gas: 
Bragg reflections of a quantum many-body wavepacket 
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The dynamics of the momentum distribution and the reduced single-particle density matrix 
(RSPDM) of a Tonks- Girardeau (TG) gas is studied in the context of Bragg-reflections of a 
many-body wavepacket. We find strong suppression of a Bragg-reflection peak for a dense TG 
wavepacket; our observation illustrates dependence of the momentum distribution on the interac- 
tions/wavefunction symmetry. The momentum distribution is calculated with a fast algorithm based 
on a formula expressing the RSPDM via a dynamically evolving single-particle basis. 
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The possibility of constraining atomic gases to one- 
dimensional (ID) geometries P, has lead to experi- 
mental realizations of exactly solvable ID models describ- 
ing interacting bose gases [J, [f|. At low temperatures, 
low linear densities, and strong repulsive effective inter- 
actions, these ID atomic gases enter a Tonks- Girardeau 
(TG) regime @, 0, 3, which is described by an exactly 
solvable model of ID bosons with "impenetrable core" re- 
pulsive interactions [H . Two recent experiments achieved 
the TG regime and observed the properties of a TG gas 
0, ||. One particularly interesting aspect of these ID 
systems is their nonequilibrium dynamics. A recent ex- 
periment studying nonequilibrium dynamics of a ID in- 
teracting bose gas (including the TG regime) has shown 
that its momentum distribution does not need to relax 
to thermodynamic equilibrium even after numerous colli- 
sions @. These experimental advances and the possibil- 
ity of exactly solving the TG model 0, [l(| , motivate us 
to study the momentum distribution of the dynamically 
evolving TG gas. 

The TG model is exactly solvable via Fermi-Bose map- 
ping, which relates the TG gas to a system of noninteract- 
ing spinless ID fermions Many properties of the 

two systems such as the the single-particle (SP) density 
0) Q3 or the thermodynamic properties [ll[ are iden- 
tical. However, quantum correlations contained within 
the reduced single-particle density matrix (RSPDM), or 
the momentum distribution of the TG gas ns(fc), con- 
siderably differ from those of the ideal Fermi gas %(fc) 

[IlIiSIiEliEIiEISIiEElIMIlE Hp. Although 

the exact many-body wavefunction describing TG gas 
can be written in compact form [H, [l(| , the calculation of 
the RSPDM and the momentum distribution is a difficult 
task [H, [3, [H, H H EI M, HH, HI . In the station- 
ary case, the RSPDM and rae(fc) were studied for a TG 
gas on the ring [l3|, [l8| and in the harmonic confinement 
flil [l5l fl7l. Il8j . In the homogeneous case, the momentum 
distribution has a singularity at k = 0, ns{k) oc k~ x l 2 
(HI, and slowly decaying tails ns{k) oc k~ 4 [lg]. In both 
the homogeneous and the harmonic case, the occupation 
of the leadingnatural orbital (effective SP state) is oc \HV 
for large N [18j |. An analytic approximation for momen- 
tum distribution of a TG gas in a box has been made by 



generalizing the Haldane's harmonic-fluid approach [l6| 
In the time-dependent case, the RSPDM and momen- 
tum distribution of the TG gas was studied in a har- 



monic potential with the time-dependent frequency 21 1 ; 
dynamics was solved with a scaling transformation 21 1 . 
Irregular motion, and the dynamics of the momentum 
distribution, was studied numerically for different inter- 
action strengths (up to the TG limit) in Ref. [19j; solu- 
tions for N — 6 bosons were presented. Several recent 
studies have addressed the dynamics of hard-core bosons 
(HCB) on the lattice [H H4]. Numerical studies of this 
model revealed a number of interesting results including 
fermionization of the momentum distribution during ID 
free expansion [20], and the possibility of relaxation of 
this system to a steady state, which carries memory of 
the initial conditions [22j. However, the behavior of the 
discrete HCB-lattice model is not equivalent to the TG 
bosons in a continuous potential 23] . A feasible numeri- 
cal study of the RSPDM and related observables during 
motion in a continuous potential V(x. t) demands an ef- 
ficient method for the calculation of the RSPDM, inde- 
pendent of the external potential, the state of the system, 
and which would be operative for a larger number of par- 
ticles. 

Here we study dynamics of the momentum distribu- 
tion, the RSPDM, natural orbitals (NOs), their occupan- 
cies, and Shannon entropy for a TG gas in a continuous 
potential. Our calculation is based on a formula express- 
ing the RSPDM via a dynamically evolving SP basis; the 
method does not depend on the external potential, the 
state of the system, and it is operative for a larger number 
of particles. The method is employed in studying Bragg 
reflections of a TG many-body wavepacket in periodic po- 
tentials. A comparison of the TG bosonic (ub) and non- 
interactiong fermionic (rip) momentum distributions il- 
lustrates the influence of interactions/wavefunction sym- 
metry on this observable. The momentum distribution 
of the ideal fermi gas displays a beating peak at the edge 
of the Brillouin zone. In contrast, such a Bragg-reflection 
peak is completely absent for a dense TG wavepacket. As 
the TG wavepacket reflects from the potential, it under- 
goes a rapid decrease of coherence, characterized by the 
increase of entropy and decrease of spatial correlations. 



2 



The model- We consider dynamics of N indistin- 
guishable bosons in ID configuration space, located 
in an external potential V(x,t), and interacting via 
impenetrable pointlike interactions [J. The bosonic 
many-body wavefunction ip B (xi, ■ ■ ■ , %n, t) describing 
the state of this system is related to a fermionic 
wavefunction ipp(xi, . . . , xjsr, t), which describes a 
system of N noninteracting spinless ID fermions: 
ip B (xi,...,x N ,t) = A(x 1 ,...,x N )ip F (xx,...,x N ,t), 
where A = ni<i<j<Arsgn(a;i — Xj) is a "unit anti- 
symmetric function" ; this is the famous Fermi-Bose 
mapping (ij. The dynamics of the fermionic wave- 
function ipp can be constructed from the Slater de- 
terminant ipp(xi, . . . , Xn, t) = y/l/Nl det[ip m (xj , t)], 
where tp m {x,t) denote N orthonormal SP wavefunctions 
ip m (x,t) obeying 



ih 



at 



ip m (x,t), m = 1, 



,N. 
(1) 

The exact many-body wavefunction of the TG system is 



where the entries of the matrix P are Pij(x, y, t) — 5ij — 
2 f dx'ijj*{x', t)t/jj(x', t); we have assumed x < y without 
loss of generality. 

Derivation of formula ([J) as follows. Define permuta- 
tions (k 2 . . . k N ) = P{1 . . . i - 1 i + 1 . . . N), {h ■ ■ ■ In) = 
Q(l . . . j — 1 j + 1 . . . N), and their signatures e(P) and 
e(Q). From the definition of the RSPDM and Eq. © it 
follows that 



A, 



N 



^y— ^ J [[dx n sgn(x - x n )sgn(y - x n ) 

J2e(P)r k2 (x 2 )...r kN M 
p 

^2 e(Q)'0i 2 (x 2 )--- ipi N (xn) 



(-1) 



(iV-l) 
(-l) i+J detP y , 



(5) 

(6) 
(7) 



ip B = A(xi, . . .,X N ) 



1 N 

det UPmixjM (2) 

iv ! m,j— 1 



i.e., its evolution is constructed after solving Eq. (TT]). 

The RSPDM of the TG system, p B (x,y,t) = 
N Jdx 2 ■ ■ ■ dx N ip B (x,x 2 ,..., x N , t)*ip B (y, x 2 ,..., x N , t), 
furnishes the expectation values of one-particle 
observables such as the position density 
p B {x,x,t), or momentum distribution n B (k,i) = 
(2Tr)- 1 Jdxdye ik ^-^p B (x,y,t) |13J. The natural or- 
bitals (j>i(x,t) (NOs) of the TG system, obtained as 
eigenfunctions of the RSPDM, 



dyp B (x,y,t)(f)i(y,t) = \i{t)(j)i{x,t), i= 1,2, 



represent effective SP states, while eigenvalues Xi(t) 
represent their occupancies [lj]. The SP wavefunc- 
tions ipm (xj , t) are NOs of the fermionic system, with 
occupancy unity, because the fermionic RSPDM is 

pp(x,y,t) = ELiC^-Wmtl/,*) 0- The momen- 
tum distributions can be expressed via the Fourier trans- 
form of the NOs, np(k,t) — 53m=i li'mik, t)| 2 , and 

The method.- The RSPDM can be expressed in terms 
of the dynamically evolving SP basis: 



PB 



(x,V,t) =^2ijj*{x,t)A lj (x,y,t)ip : j(y,t). (3) 



The N x N matrix A(a;, y, t) = {Aij(x, y, t)} is 



A(x,y,t) = detP(P" 



-1\T 



where Py is a minor of matrix P obtained by crossing 
its ith row and jth column. Equation ([6|) is obtained 
after rearranging the product factors of Eq. ([5]), and 
formally performing the integrations Pk n ,i„ = <5fc„,;„ — 
2 \l dx'^* kn (x', t)i/j ln 0', t). Eq. follows from the 'defi- 
nition of a determinant [18|] ■ Eq. ^ follows immediately 
from Eq. ([7]) and the formula for the matrix inverse via 
algebraic co-factors. 

Dynamics in the periodic potential- The richness of 
the dynamics of ultracold bose gases in optical lattices 
[24j motivate us to numerically (exactly) study the evo- 
lution of a quantum many-body wavepacket in a con- 
tinuous periodic potential V„[x) = V p (x + D) (also re- 
ferred to as the lattice); periodic boundary conditions 
are assumed, i.e., dynamics occurs on a ring of length 
L = n s D. The gas (wavepacket) is initially localized 
within a region significantly smaller than L, and it is 
given a certain amount of momentum. During dynam- 
ics, the many-body wavepacket will disperse on the ring. 
The dynamics of the TG momentum distribution n B is 
affected by the exchange of the momentum between the 
lattice and the gas, the many-body interactions, and the 
bosonic symmetry of the wavefunction. On the other 
side, the related fermionic momentum distribution n F is 
affected by the lattice and the Pauli exclusion principle. 
We find it illustrative to compare time-evolution of the 
two momentum distributions, as it illustrates the influ- 
ence of the interactions/wavefunction symmetry on this 
observable. 

In our numerical simulations we consider motion of 
87 Rb atoms in the potential V p (x) = Vo cos 2 (irx / D) , 
where D = 391.5 /mi, and Vo — 11.9 eV unless spec- 
ified otherwise; n s — 52. For concreteness, we assume 
that the SP wavefunctions describing the wavepacket at 



(4) t = are ip m (x, 0) 



i(x)e 



ih' X 



, N. where 
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FIG. 1: (color online) Dynamics of momentum distributions, 
(a) riF(k,t) in the initial stage of the evolution (down), and 
after long time propagation (up); beating at k = —tt/D is 
a signature of Bragg reflections, (b) Initial momentum dis- 
tributions riF(k, 0) and n_e(fc,0) for N = 25 bosons (solid 
lines); dot-dashed line depicts X^ m =i IVK^i *)| 2 ( u p)> an d 
5Zm=i ^Wl'/K^' *)| 2 (down) at f = 34.5 ms (the area beneath 
the dot-dashed curves up and down is equal), (c) ns(fc,t) for 
N = 25 bosons in the initial evolution stage (down) , and after 
long time propagation (up) ; the signature of Bragg reflections 
at k — —tt/D is absent, (d) Dynamics of ns(fc, i) for N = 3 
bosons; there is beating at k = —tt/D. 



u m is the mth SP eigenstate of the harmonic poten- 
tial Vh(x) = mw 2 /2, u> = 2tt 316Hz; such a many-body 
wavepacket corresponds to a ground state of the gas in 
harmonic confinement, with a momentum k' per parti- 
cle imparted to the wavepacket. The initial expectation 
value of the SP momentum k' = JknB(k 7 0)dk is chosen 
to be exactly at the edge of the Brillouin zone k' — tt/D. 
Although such an excitation is non-trivial to prepare, 
current high level of experimental techniques [l|, 0, M 0] 
strongly suggests that it is more than just a theoretical 
curiosity. 

It should be emphasized that the expectation value 
of the SP momentum is identical (at all times) for 
TG bosons and noninteracting fermions, = 
Jdkkns(k,t) = fdkknp(k,t) — (k)p. Nevertheless, 
their momentum distributions show remarkable differ- 
ences. Figure Q3a) shows rip in the initial stage of the 
evolution, and after long-time propagation (when the gas 
is already well-dispersed over the ring). A sharp peak 
beating up-down at the edge of the 1st Brillouin zone 
k = —tt/D arises from Bragg reflections. The fermionic 
momentum distribution is np{k,t) — J2 m =i ; 
a few of the SP spectra \ip m (k,t)\ 2 are initially overlap- 
ping the edge of the Brillouin zone at tt/D; as the dy- 
namics of ipmikjt) are uncoupled, the spectra \ip m (k, t)\ 2 
of those NOs display a beating Bragg-reflection peak at 
—tt/D [see dot-dashed curve in Fig. [TJb) up], which is 
reflected onto np(k, t). 




FIG. 2: (color online) Dynamics of bosonic NO occupations 
Xi and entropy S. (a) Xi(t) at times t = 0, 0.40, 2.15, 33, 34.5 
ms; lattice depth is Vb = 11.9 eV. (b) The entropy S(t) 
for three different lattice depths Vb = 0, 4.77, 23.85 eV. (c) 
\Pf{0, x,t)\ and \pb(0, x,t)\ at t = (solid line) and t = 34.5 
(dot-dashed line); Vb = 11.9 eV. (d) The Bloch-wavevector 
distribution for 25 bosons at t = 0, and t — 34.5 ms; upper 
picture (lower picture) corresponds to lattice depth Vb = 2.36 
eV (Vb = 11.9 eV), respectively. See text for details. 



The bosonic momentum distribution tib at t = is 
shown in Fig. HJb). The sharp peak of ns(fc, 0) is located 
exactly at the edge of the Brillouin zone, and it is much 
sharper than the peak of np{k 1 Q). From this one may 
erroneously conclude that there would be a sharp beating 
peak originating from Bragg reflections at —tt/D. How- 
ever, this signature of Bragg reflections is absent. This 
is illustrated in Fig. [IJc), which shows a contour plot of 
the bosonic momentum distribution ns{k, t) for N = 25 
bosons propagating in the potential V p . The signature 
is absent both at the beginning of the motion, when the 
wavepacket is still localized, and after it spreads over the 
ring. In the long tome propagation tib collectively os- 
cillates due to the momentum-exchange with the lattice 
i(k)B = (^)f)j but the changes in its shape are small. 
Our simulation clearly depicts that when the momen- 
tum is being transferred by the lattice to the TG gas, it 
redistributes among bosons; this leads to a smooth distri- 
bution without a beating Bragg-reflection peak. Unlike 
the fermionic NOs, the low-order bosonic NOs do not 
display Bragg-reflection peaks due to strong (nonlinear) 
coupling arising from interactions [see dot-dashed curve 
in Fig. [TJb) down]. 

However, Bragg-reflection peaks can be obtained for a 
smaller density of the TG gas. This is illustrated in Fig- 
ure [ljd) showing an identical numerical simulation but 
with N = 3 bosons. In this CclS6, cLS bosons disperse 
on the ring, their density is sufficiently low, leading to 
the fermionization of the bosonic momentum distribu- 
tion discussed in Ref. [2(J; consequently one observes a 
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beating peak at k = —ir/D. 

Further insight into Bragg reflections of the TG many- 
body wavepacket follows from the behavior of the bosonic 
NOs (f>i and their occupancies A^. Initially (t = 0), 
a few of the leading NO occupancies are fairly large 
[see Fig. [2(a)], which is characteristic for a cold Bose 
gas. However, when the evolution begins, the low or- 
der XiS rapidly decrease, while the number of NOs with 
non-negligible occupations increase [Fig. [2(a)] . Fig- 
ure (Kb) illustrates the time-evolution of the Shannon 
entropy S(t) = -EiPi lo gP» where Pi(t) = \ l (t)/N, 
for different lattice depths Vo; the entropy S increases 
faster, and saturates at a higher value for a deeper lattice. 
Figure EJc) shows bosonic (fcrmionic) quantum correla- 
tions fx B (x,x',t) = pb{x, x' , t)j p B {x, x, t)p B (x', x', t), 
[Hf{x, x' , t), respectively] at t = and t = 34.5 ms; 
in contrast to p B , the correlations of nonintcracting 
fermions are not considerably changed during evolution. 
Figures [2(a)- (c) clearly illustrate the dynamical loss of 
coherence of the TG wavepacket, which is more rapid 
for a deeper lattice. This results from the interplay of 
the many-body interactions and scattering from the lat- 
tice. Namely, interactions couple bosonic NOs thereby 
providing a mechanism for the time-change of their oc- 
cupancies. For a deeper lattice, the initial wavepacket ef- 
fectively excites a larger number of system's eigenstates; 
this is illustrated in Fig. [2(d) which shows the diago- 
nal of the RSPDM represented in the Bloch-wave basis 



(extended Brillouin-zone scheme) for two different lattice 
depths. For a deeper lattice, the dynamics effectively in- 
volves a larger number of frequencies (i.e., energies), and 
it is more irregular. 

Before closing, we note that the dynamics of the TG 
gas is related to the paraxial propagation of partially- 
incoherent light (PIL) beams in linear ID photonic 
structures [251 ] . Furthermore, the behavior of partially- 
condensed weakly-interacting Bose gases is similar to PIL 
in noninstantaneous nonlinear media [26| . These facts 
motivate us to explore the recently observed phenomena 
with incoherent light in photonic lattices 27J , within the 
context of quantum-dynamics of interacting bosons. 

In conclusion, we have studied dynamics of the mo- 
mentum distribution, RSPDM correlations, natural or- 
bitals and their occupancies, and the entropy of the TG 
gas out of equilibrium. We analyzed Bragg reflections 
of the TG many-body wavepacket and found that their 
signature (observed as a beating resonant peak in the 
momentum distribution of the corresponding noninter- 
acting fermionic gas) may be considerably suppressed by 
the TG many-body interactions. We have employed a 
fast numerical method, applicable for versatile continu- 
ous potentials, and operative for larger number of parti- 
cles. Our results open the way for further studies of the 
RSPDM and related observables of the TG gas, both in 
the static and time-dependent cases. 
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